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The extended gaussian ensemble (EGE) is introduced as a generalization of the canonical en- 
semble. The new ensemble is a further extension of the Gaussian ensemble introduced by J. H. 
Hetherington [J. Low Temp. Phys. 66, 145 (1987)]. The statistical mechanical formalism is de- 
rived both from the analysis of the system attached to a finite reservoir and from the Maximum 
Statistical Entropy Principle. The probability of each microstate depends on two parameters (3 and 
7 which allow to fix, independently, the mean energy of the system and the energy fluctuations 
respectively. We establish the Legendre transform structure for the generalized thermodynamic 
potential and propose a stability criterion. We also compare the EGE probability distribution with 
the q-exponential distribution. As an example, an application to a system with few independent 
spins is presented. 
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I. INTRODUCTION 



The development of statistical mechanics based on en- 
semble theory is founded on the postulate of "equal a 
priori probabilities" , which is assumed to apply to all 
microstates consistent with the given macrostate of an 
isolated system pj . The corresponding statistical ensem- 
ble is the so-called microcanonical ensemble. A repre- 
sentative system in this ensemble has all "mechanical" 
variables such as energy E, volume V, magnetization M 
etc., fixed. For convenience in calculations, other ensem- 
bles are used which invariably suppose the existence of 
a subsidiary system or reservoir in contact with the ac- 
tual system. For instance, in the canonical ensemble the 
walls of the system permit an exchange of energy with 
the reservoir while in the grand canonical ensemble, both 
energy and matter can be exchanged. In general, the dif- 
ferent ensembles are constructed by allowing one or more 
mechanical variables to fluctuate. The exchange of each 
of these variables is controlled by a parameter which is a 
characteristic of the reservoir. For instance, in the case 
of the canonical ensemble, this parameter is precisely the 
temperature of the reservoir and determines the mean 
energy of the system. Actually, this is adequate when 
the reservoir is a very large system that can exchange 
arbitrary amounts of energy, without modification of its 
intensive properties. In practical situations, this is not 
always the case. However, very few studies have been de- 
voted to analyse the consequences of possible deviations 
from these ideal reservoir properties. 

In this paper, we develop the statistical mechanics of 
a system that can exchange energy with a finite reservoir 
characterized by two parameters: (3 and 7. These pa- 
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rameters control independently the mean energy of the 
system and its energy fluctuations, respectively. The 
corresponding statistical ensemble represents a general- 
ization of the canonical ensemble and will be called the 
extended gaussian ensemble (EGE). A similar ensemble 
was already developed, in a more restricted framework, 
by Hetherington 2] . The author considered that the sam- 
ple system was in contact with a finite reservoir with size 
dependent properties. The so-called gaussian ensemble 
was introduced so that it is equivalent to the canonical 
ensemble in the limit of large systems, except in the en- 
ergy range of a first-order transition. Interestingly, it en- 
ables a smooth interpolation between the microcanonical 
and the canonical ensembles. Taking into account these 
features, Challa and Hetherington |J0| showed the inter- 
est of this ensemble for Monte Carlo simulation studies 
of phase transitions. They demonstrated a significant 
reduction in computer time (compared to standard sim- 
ulations in the canonical ensemble) and, its adequacy for 
distinguishing second-order from first-order transitions. 
Compared to the EGE introduced in the present paper, 
the main difference arises from the fact that in the gaus- 
sian ensemble the sample and the reservoir are assumed 
to be statistically independent which implies the additiv- 
ity of the corresponding entropies. This is not assumed 
in our formalism. The consequences are important and 
will be discussed in depth in this work. 

The present formalism can be considered as an alter- 
native to the statistical mechanics based on non-additive 
generalized entropies. Actually the study of such gener- 
alized entropies has generated a lot of interest in the past 
fiveteen years. The motivation for the so called Tsallis 
statistical mechanics has been to extend the standard 
Boltzmann-Gibbs framework to include non-extensive 
systems 5J. Among different interpretations, it has been 
suggested that Tsallis formalism corresponds to an en- 
semble describing a system attached to a finite reservoir 
0, Q- Although a large number of papers have been 
published, the physical meaning of many related issues is 
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still open to discussion 0, . The EGE formalism that 
we propose in this paper, provides a clear and consis- 
tent framework for the statistical mechanics with a finite 
reservoir. 

The paper is organized as follows: in sections [H] and 
IIIII the EGE is founded from the analysis of a contact 
with a finite reservoir and from the Maximum Statistical 
Entropy Principle respectively. In section IIVI the main 
thermodynamic relations are derived. In section we 
highlight the non-additive nature of the thermodynamic 
formalism. In section lvTl a stability criterion is proposed. 
In section IVlII the equilibrium distributions of the EGE 
are compared with the g-exponential distributions. In 



and 



section [Villi we present an example of application to a 
system of independent spins. Finally, in section IIXI we 
summarize and conclude. 



II. CONTACT WITH A FINITE RESERVOIR 



d 2 S 2 



dEi 
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The standard canonical ensemble is characterized by a 
reservoir with constant f3 (independent of E2), which im- 
plies 7 = and there is no term beyond the first order 
term in the Eq. However, in the present paper, we 
consider a more general reservoir for which 7^0. Thus 
in this EGE, the reservoir is characterized by the pair 
of parameters (3 and 7. The thermodynamic meaning 
of these parameters will be clarified in the following sec- 
tions. To explicitly highlight the effects of this modifica- 
tion and also for the sake of simplicity, we assume that 
the cubic and higher order terms vanish. Then substitut- 
ing Q in © and denoting the energies of the microstates 
of the sample by e, [i = 1, M), we obtain 



Let us consider a system (that we will call the sample) 
in contact with a reservoir. Let us call the energy of the 
sample E\ and the energy of the reservoir E 2 . The sample 
and the reservoir together form an isolated system so that 
E = Ei +E2 is constant. Let us also define il 2 (E 2 ) as the 
number of microstates of the reservoir. Following Callen 
[To), the probability that the system 1 is in a certain 
microstate with energy E\ is given by 



n 2 (£-£i) 
n 1+2 (E) ' 



(i) 



where Q,i +2 (E) is the total number of states available for 
the set 1 + 2 (Note that we do not assume that Oi + 2 
can be factorized as a product U,\^l 2 )■ Let us define the 
entropy of the reservoir as: 



S 2 {E 2 ) = In Sl 2 (E 2 ), 



(2) 



(all along the paper we choose the entropy units so that 
Boltzmann's constant, ks = 1.) Therefore: 



e S 2 (E-E 1 ) 

n 1+2 (E) ■ 



(3) 



The energy of the sample will, in general fluctuate. Let 
us call U its mean (equilibrium) value. We can develop 
S 2 (E — Ei) around the equilibrium value E — U as: 



S 2 {E-E X ) = S 2 {E-U) + 



1 d 2 S 2 



dth 
dE 2 



{U-E x ) 
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(4) 



The derivatives in the right hand side of this expression 
are quantities which depend only on the reservoir. We 
define 



dS 2 
dE 2 



= 0, 



(5) 



Pi = — exp[-/5e l - 7 (e l -C/) 2 ], (7) 

AG 

where the normalization constant Zq is given by 



M 



(8) 



The subscript G only indicates the "gaussian" form of 
the probabilities. Note that U is the mean energy and 
must be obtained self-consistently from the equation: 



M 

UZ G = y^exp[ 



-pet- 7 (e* - Uf 



(9) 



Eqs. 0, © and iJjJJ reduce to the standard canonical 
ensemble definitions when 7 = 0. Therefore, it is natu- 
ral to relate the parameter 7 with the finite size of the 
reservoir. 



III. MAXIMUM STATISTICAL ENTROPY 
PRINCIPLE 

In this section, we derive the probability law of Eq. Q 
from different arguments. This leads to a better under- 
standing of the parameters /3 and 7 as parameters char- 
acterizing the "equilibrium" distribution of the sample. 
To derive the probability distribution from the Maximum 
Statistical Entropy Principle, we maximize the standard 
Gibbs-Boltzmann-Shannon entropy given by 



M 



Sg = ~y) J?, In 



(10) 



E-U 



subject to the constraints of normalisation of the proba- 
bility, the given mean value of the energy and the fixed 
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value of the fluctuations, respectively as 

M 



E* = ( n ) 

M 

»=i 

(( ei -C/) 2 )=^( ei -«7) 2 Pi = W. (13) 



Then the maximization procedure is done by introducing 
the Lagrange multipiers A, (3 and 7 for the respective 
constraints, and maximizing the following functional C: 

c = -Eft - a ^Eft - M - z 3 (e eift - ^ 
-7 ^E^-^'ft-^) • 



By requiring the condition: 



OPi 



(14) 



(15) 



it is easy to see that the optimum form of the proba- 
bility distribution is given by the expression in Eq. J7J. 
Therefore f3 and 7, within this context, are simply La- 
grange multipliers that allow to fix, self-consistently, a 
mean value of the energy U — (e,) and a specific value of 
the variance W — ((ej — U) 2 ). 



IV. THERMODYNAMIC RELATIONS 

We define a thermodynamic potential $(/3, 7) as 

$00,7)= In Z G . (16) 

By differentiating Eq. it can be straightforwardly 
obtained that: 



d$\ 



The second derivative renders: 
'<9 2 $\ fdU 



8PJi W/ 7 M^- 1 (/3,7)-2 7 



(17) 
(18) 

. (19) 



which represents a generalization of the standard formula 
for energy fluctuations in the canonical ensemble. It is 
natural to define the extended heat capacity as: 



a i' du \ P 2w 



d{i) 1 1-27W 



(20) 



This equation is the same that was already derived in 
Ref.0. Note that, contrary to what happens in the stan- 
dard canonical ensemble, the positivity of the fluctua- 
tions W does not guarantee the positivity of C. 

For 7 — > 0, it is seen that the relations (|17fl and l|19f> go 
to the corresponding relations for the case of canonical 
ensemble. Also in this limit, from Eqs. (|18fl and l|19|) we 
get an interesting relation given by: 

' d 2 $ s 

w 



(21) 



which resembles in form with a diffusion equation. 

The entropy Sg as given by (Unjl is the inverse Legendre 
transform of $(/?, 7), and can be expressed as: 

S G (U,W) =(3U + 7 iy + $, (22) 

whereby Sq is a function of the specified values of the 
constraints i.e. U and W. Therefore we have the follow- 
ing thermodynamic relations 

'dS G 



9U /W 
fdS G \ 

\dw) u 



= (3, 



(23) 
(24) 



V. NON-ADDITIVITY 

We remark that although the thermodynamics of a sys- 
tem in the EGE is well defined by the equations in the 
previous section, it is not straighforward to establish a 
mutual equilibrium condition for two different systems 
that would allow to establish a zeroth law (or, equiv- 
alently, an intensive temperature) [llj . This problem is 
due to the non additive character of the potential $(/3, 7). 
Let us consider two systems 1 and 2 with hamiltonians 
Hi and H^. By applying the rules derived in the previ- 
ous sections independently to the two systems, one can 
derive the thermodynamic potentials $i(/3, 7), $2(/?, 7) 
as well as the mean energies U±((3, 7) and U2(P, 7). One 
can then try to solve the composite system 1 + 2 with 
hamiltonian Hi + Hi- It is easy to verify that the new 
potential $1+2 (/?, 7) satisfies: 

$1+2(0,7) = $i(/3,7) + $2(0,7) 



where L7+2 is the mean energy of the composite sys- 
tem. Note that U1+2 as well as the average indicated 
by the angular brackets are computed with the probabil- 
ity distribution corresponding to the composite system 
1 + 2 which, in general, cannot be written as a product 
of probability distributions for systems 1 and 2. 

The average values are, in general, non-additive 
(U1+2 i= Ui + U2) ■ But even if additivity of U is im- 
posed, the potential $ remains non-additive. The cor- 
rection term depends on the microscopic details of the 
two hamiltonians Hi and H2- 
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This lack of additivity does not allow to define an 
equivalence relation of "mutual" equilibrium. Consider 
that two systems are known to be independently in equi- 
librium with a bath characterized by (3 and 7. On putting 
the two systems in contact with bath simultaneously, 
their properties (probability distribution) will change. 
The parameters and 7 are, consequently, properties 
of the bath and cannot be considered as properties of the 
system, contrary to what occurs with the corresponding 
/3 in the standard canonical ensemble. 



VI. STABILITY CRITERION 

In standard thermodynamics, the stability criterion 
— dU /d[3 > is derived from the condition of maximum 
entropy. The derivation |lfl | considers a partition of an 
isolated system into any two subsystems. By allowing 
the two subsystems to alter their energies at fixed to- 
tal energy, one can analyze the entropy change when the 
system is (virtually) displaced out of equilibrium. The 
condition of the maximum total entropy allows to de- 
duce that the equilibrium state corresponds to a state 
with homogeneous (3 (equilibrium condition) and with 
—dU/d/3 > (stability condition). For the derivation, 
nevertheless, additivity of the entropy of the two subsys- 
tems must be used. 

Within our new formalism, an additivity assumption 
cannot be used. Therefore it is not straightforward to 
establish a stability criterion. Although we cannot give a 
rigorous proof, in this section we provide some evidences 
that the same criterion (—dU/d(3 > 0) must hold. 

First of all, it interesting to remark that the require- 
ment of Eq. © that allows to find U can be rewritten 
as an extremal condition. Consider the definition of an 
"out-of-equilibrium" potential: 



9(p,>y,U)=hxZ. 



(26) 



Note that here we are considering (3, 7 and U as indepen- 
dent variables. U shall not be regarded as the internal 
energy but as a parameter that allows virtual displace- 
ments out of equilibrium. It can be checked that the self 
consistent Eq. © can then be written as: 



0. 



(27) 



Therefore, only the extrema of W correspond to equilib- 
rium solutions: U = U(f3,~f). By substituting in \I> one 
gets the equilibrium potential: 



$(/3, 7) =*(/?, 7,tf(/?,7)). 



(28) 
±00, 



Secondly, from Eq. © note that if 7 > and U 
then Zq — > and therefore W — > —00. 

We can use this "out-of-equilibrium" potential to de- 
fine a stability criterion. It is straighforward to compute 



its second derivative: 



a 2 * 



2 7 {2jW - 1) . 



(29) 



/3>7 



Note that the positivity of C in Eq. (121 H would guarantee 
that this second derivative of \& is negative and there- 
fore the state of equilibrium corresponds to maxima of 
^(0, 7, U) with respect to U displacements. Although, 
contrary to what happens in the standard canonical en- 
semble, the positivity of W does not ensure the sign of 
—dU/d/3 in general, at least we can derive that for small 
and positive values of 7: 



< - 



dU 
dp 



d 2 ^ 
dU 2 



< 0. 



(30) 



0,1 



Thus we suggest that this is the stability criterion to be 
used within the EGE, and we will use it in Section lYlIII 
for the analysis of some examples. 



VII. COMPARISON WITH g-EXPONENTIAL 
DISTRIBUTIONS 

^-exponential distributions are the central predictions 
of the generalized statistical mechanics proposed by Tsal- 
lis |5| . These distributions have been considered as model 
distributions to describe various complex systems at their 
stationary states 

El 13 El- The general form of such 
distributions is given by p(x) ~ e q (x), where the q- 
exponential is defined as e q (x) = [1 + (1 — q)x] 1 ^ 1 ^ q ^ . 
This function goes to the usual exp(x) function for q — > 1. 
For defmitcness, we restrict to the range < q < 1. 

In this section, we compare the ^-exponential distri- 
butions with the equilibrium distributions of the EGE. 
But first, we show how to derive the g-exponential distri- 
butions by generalizing the canonical ensemble approach, 
along the lines of section^] We define a parameter which 
is in general, a function of the energy E2 of the reservoir 



(31) 



At equilibrium, it attains the value given by Eq. (|5jl. We 
impose that (3(E 2 ) satisfies: 



d ( 1 



dE 2 \/3(E 2 ) 



Q, 



(32) 



where Q is a positive valued constant. From Eqs. i|32|) 
and H31fl . we obtain 



d R( TP \ - d2S2 - '1 <- 



QP 2 {E 2 ) 



In general, for all integer values of n 



dE. 



± = {n-l)\{-QT-'(3 n {E 2 ). 



(33) 



(34) 
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Now unlike in Eq. if in the expansion of S(E — E\) 
around the equilibrium value [E — U), we retain deriva- 
tives of S 2 upto all orders, then we have 



°° -i Ana 

S 2 (E-E 1 )=S 2 (E-U) + Y / - i " 



n! dES 

n=l 2 



(C/-^i) n . 
(35) 

On applying Eq. I|34|) for the case of equilibrium, we can 
write 

oo 

s 2 (e - e 1 ) = s 2 (£ - u) + -{-QT^^iu - E x ) n , 



n=l 



(36) 

where note that (3 is given by its value at equilibrium. 
The equilibrium probability distribution is then given 
from © as 



p(Ei) ~ exp 



00 i 

v-(-Qr-\3 n (f/-£i) r 

^ — ' rj 



.n=l 



(37) 



To compare Eq. (|37|l with the q-exponential distribution 
given by 



P g (Ex 



\J3{U-E X )] = [1 + (1 - q)0(U - Ei)} 



1/(1-9) 

(38) 



(39) 



we rewrite the g-exponential as 

\RtTT j? W [ Ml + (l-g)^-^i)l 
e q \fi{U-Ex)\ =exp _ ^ 

and expand the In function using the series ln[l + x] = 

2 3 4 

x-^T + ^j — + provided that — 1 < x < 1. Thus 
we can write 

e q [P(U-Ei)} = exp 



oo 1 

^-{-(i-e)}"- 1 /? 71 ^-^; 

^ — ' rj 



.n=l 



(40) 

for -1 < (1 - q)j3(U - E t ) < 1. Thereby, on identify- 
ing Q = (1 — (?) we may say that the general equilibrium 
distribution of Eq. I|37|l based on assumptions (|31|l and 
(|32|l . is identical to a g-exponential distribution. Assum- 
ing that the relevant q values are quite close to unity, 
we may keep terms only upto second order as done in 
Eq. Then the equilibrium g-distribution for system 
1 being in microstate i of energy ej can be written as 



-Pe i --(l-q)l3 2 (e i -U) 2 



(41) 



where Z q is the normalization constant. 

On the other hand, for the case of EGE, instead of fix- 
ing the derivative of (3~ x (Eq. I)32[0. we fix the derivative 
of P as follows: 



dE- 



■(3{E 2 



-2 7 , 



(42) 



where 7 is independent of E 2 . This ensures that the 
higher order [n > 2) derivatives of S 2 vanish. On com- 
paring (|41|l and (7J, we note that (1 — q) plays the role 
analogous to 7. 



It may be remarked that if we identify parameter 
P{E 2 ) = 1/T(E 2 ) as the inverse temperature, then 
Eq. 1)32(1 implies that the heat capacity of the reservoir 
C 2 = dE 2 /dT = Recently, the ^-exponential distri- 

butions have been discussed in the context of a reservoir 
with finite heat capacity [3 ■ On the other hand, following 
Gibbs' approach to the canonical ensemble, but instead 
using the ^-generalized Boltzmann entropy, ^-exponential 
distributions were derived in Ref. 



VIII. APPLICATION TO A SYSTEM OF 
INDEPENDENT SPINS 

A. Single Spin 

As a first example of the EGE, we apply our formalism 
to the problem of a system with only two energy levels. 
Let us consider a single spin s = ±1 in the presence of a 
constant external magnetic field B. The hamiltonian of 
the system reads: 

H = -Bs. (43) 
The partition function is given by: 

Z G = e PB e -^-B-Uf + e - P B e - l{ B-Uf t (44) 

where the mean energy U is the solution of the self- 
consistent equation: 

U = -Be' m e-^- B - u f + Be-V B e-^-V? . (45) 

The dependence on B can be easily overcome by defining 
the reduced units: 

U* = U/B, (3* = f3B, 7* = 7B 2 . (46) 

Thus, Eq. lj4"5)l becomes: 

(47) 



The numerical solution of this equation is plotted in 
Fig. ^ The behaviour of U* as a function of 1/(3* is 
shown for different values of 7*. For 7* = 0, one recov- 
ers the behaviour U* = tanh(/3*) corresponding to the 
case of a system in contact with a infinite reservoir. For 
7* ^ 0, U* is smaller, indicating that it is more dificult to 
disorder the system by decreasing (3* . It is interesting to 
note that for 7* ~ 0.5, there is a change in the behaviour 
at (3* — > 0. Above this value of 7*, the system is not 
able to disorder completely anymore and always keeps a 
certain magnetization (m = (s) = —U*). This can be 
regarded as a "phase transition", that occurs at (3* = 0. 
This change in the behaviour occuring at 7* ~ 0.5 can 
also be seen by plotting the entropy S as a function of 
1/(3* for different values of 7*. This is shown in Fig. 
For 7* < 0.5 the entropy tends to In 2 for (3* — > 00, 
whereas it tends to a lower value for 7* > 0.5. 
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FIG. 1: Behaviour of the mean reduced energy U* as a func- 
tion of 1/(3* for several values of 7* in a system of a single 
spin. 
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FIG. 2: Behaviour of the entropy S as a function of 1/(5* for 
several values of 7* in a system of a single spin. 



B. Two spins 

As a second step, it is also very instructive to study a 
system of two independent spins. This will illustrate the 
non-extensive behaviour of the solution. In this case, the 
numerical solution of the self consistent equation for 
the mean energy renders the behaviour shown in Fig. [3] 
For the values of (3* and 7* for which more than one so- 
lutions are possible, we have used the stability criterion 
proposed in Section IVII to decide which is the "equili- 
birum" solution. As can be seen, for 7* > 0.49 a discon- 
tinuity occurs associated with a sudden loss of order in 
the system. Although the system is far from the thermo- 
dynamic limit, this change shares many similarities with 
a phase transition. Fig. 0] displays the behaviour of the 



corresponding energy fluctuations. It can be seen that 
W* = ((H 2 ) - U 2 )/B 2 exhibits a cusp at the transition 
for 7* ~ 0.49. For larger values of 7*, the fluctuations ex- 
hibit a discontinuity. The discontinuities are associated 
with first-order phase transitions that display metastable 
behaviour. As an example, in Fig. we show the de- 
tailed behaviour of U* as a function of for 7* = 0.6. 
In the range 2.42 < 1/(3* < 2.86, the numerical analysis 
of the self-consistent equation renders three solutions. 
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FIG. 3: Behaviour of the mean energy U* as a function of 
1/(3* for several values of 7* in a system of two independent 
spins. 
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FIG. 4: Behaviour of the energy fluctuations as a function of 
1/(3* for several values of 7* in a system of two independent 
spins. 

By analyzing the behaviour of the potential 
^((3*, 7*, U*), shown in Fig. it is easy to verify 
that two of such solutions are stable (correspond to local 
maxima of VP) whereas one is unstable (corresponds to 
a local minimum of VE" and is not plotted in Fig. [SJ. 
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FIG. 5: Behaviour of U* as a function of 1/(3* for a system 
of two spins with 7* = 0.6. The continuous line represents 
the equilibrium solution (with maximum ^) and the dashed 
lines represent metastable solutions which correspond to local 
maxima. 



The equilibrium transition jump at 1/13* ~ 2.832 is 
determined by the equality of the two maxima of 
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FIG. 6: Behaviour of the potential ty((3* , 7*, u*) for different 
values of 1/(3* for a system of two spins with 7* = 0.6. The 
equilibrium value of U* corresponds to the maximun of the 
potential <]/. 



FIG. 7: Phase diagram for a system of two spins. The con- 
tinuous line represents first-order transitions whereas discon- 
tinuous lines indicate the metastability limits. 



C. Several spins 

We have also performed a numerical study of systems 
with larger number of independent spins in the presence 
of an external field. An example is shown in Fig. [SJ cor- 
responding to a system with 4 spins (16 energy levels). 
A sequence of two consecutive phase transitions can be 
observed. As an interesting remark we want to note that 
in the case of N "non-interacting" spins s&(fc = 1, . . . N) 
in the presence of an external field, long-range forces will 
appear due to the finite size of the bath. This can be eas- 
ily seen by writing the probabilities pi for the microstates 
(i = 1, ■ ■ ■ ,2 N ) of such a system: 



p, = exp 



P 



N 

fe=i 



t N 

E 

\k=l 



s k - U* 



(48) 



Note that the development of the squared term in the 
exponent will lead to terms — ^* s^Sj wich correspond to 
antifcrromagnetic interactions among all spin pairs. A 
more detailed study of these examples is out of the scope 
of this paper. 



IX. SUMMARY 



For the system of two spins, therefore, we can plot a 
(3* — 7* phase diagram, shown in Fig. [7] The line of 
first order phase transitions ends in a "critical" point at 
(3* ~ 0.353 and 7* ~ 0.49. This point is characterized 
by the condition 1/W* = 2j* and thus, according to 
Eq. i|19fl . corresponds to a divergence of C but not to a 
divergence of the fluctuations W* , that can never diverge 
for such a system with a finite number of bounded energy 
levels. 



We have presented the extended gaussian ensemble 
(EGE) as a generalization of the standard canonical en- 
semble. The ensemble statistics have been derived by 
two methods: first by considering a system in contact 
with a finite bath, secondly from the maximum statis- 
tical entropy principle by fixing the knowledge of both 
the mean energy and the energy fluctuations. The ob- 
tained probability law depends on two parameters (3 and 
7 which are properties of the bath. Thermodynamic re- 
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rion has been suggested. Nevertheless this point as well 
as the possibility for establishing a mutual equilibrium 
criterion will need further analysis in future works. We 
have also presented an application of the EGE formalism 
to the analysis of a system of one spin and two indepen- 
dent spins. Among other interesting results, the most 
remarkable one is the possibility for occurence of a crit- 
ical point or first-order phase transitions induced by the 
finite size of the reservoir. Further, comparisons of this 
new ensemble formalism with alternative extensions of 
statistical mechanics proposed for the study of nanosys- 
tems 0, ^| or for other non-extensive systems, are 
interesting problems for research. 
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